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Abstract. We propose a novel method for calculating resonances in three- 
body Coulombic systems. The method is based on the solution of the set of 
Faddeev and Lippmann-Schwinger integral equations, which are designed for 
solving the three-body Coulomb problem. The resonances of the three-body 
system are defined as the complex-energy solutions of the homogeneous Fad- 
deev integral equations. We show how the kernels of the integral equations 
should be continued analytically in order that we get resonances. As a numer- 
ical illustration a toy model for the three-a system is solved. 

1 Introduction 

For three-body systems the Faddeev equations are the fundamental equations. 
After one iteration they possess connected kernels, consequently they are ef- 
fectively Fredholm integral equations of second kind. Thus the Frcdholm al- 
ternative applies: at certain energy either the homogeneous or inhomogeneous 
equations have solutions. Three-body bound states correspond to the solutions 
of the homogeneous Faddeev equations at real energies, resonances, as usual in 
quantum mechanics, are related to their complex-energy solutions. 

The situation becomes more complicated if the Coulomb potential enters 
into the game. The Faddeev equations were derived for short range interac- 
tions and if we simply plug-in a Coulomb-like potential they become singular. 
The solution has been formulated in the Faddeev-Merkuriev theory Q on a 
mathematically sound and elegant way via integral equations with connected 
(compact) kernels and configuration space differential equations with asymp- 
totic boundary conditions. 



*This work is dedicated to the 60th birthday of Prof. W. Glockle. 
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Recently, along the concept of a "three-potential" picture, a novel method 
has been proposed for treating the three-body Coulomb problem via solv- 
ing the set of Faddeev-Noble and Lippmann-Schwinger integral equations in 
Coulomb-Sturmian-space representation. The method was elaborated first for 
bound-state problems Q with repulsive Coulomb plus nuclear potential, then 
it was extended for calculating p — d scattering at energies below the breakup 
threshold jsj. In these calculations an excellent agreements with the results of 
other well established methods were found and the efficiency and the accuracy 
of the method were demonstrated. Also atomic bound-state problems with at- 
tractive Coulomb interactions were considered jj]. The mathematical validity 
of the procedure, in the context of the Faddeev-Merkuriev theory, has also been 
established 

In this article we seek to extend this novel method for calculating resonant 
states, the complex-energy solutions of the homogeneous Faddeev-type integral 
equations. First, along the method of Ref. ||], we recapitulate the solution of the 
homogeneous integral equations, then we show how to continue analytically the 
Green's operators onto the unphysical sheet. As an illustration of the feasibility 
of this method we calculate a model three-a system interacting via short-range 
plus repulsive Coulomb interaction. 

2 Faddeev-Merkuriev integral equations 

The Hamiltonian of a three-body Coulombic system reads 

H = H" + W + Vc,+V0 + v^, (1) 

where is the three-body kinetic energy operator, W stands for the three- 
body potential and denotes the Coulomb-like interaction in the subsystem 
a. We use throughout the usual configuration-space Jacobi coordinates Xa and 
Ua. Thus Va only depends on Xa {va — Vaixa)), while W depends on both Xa 
and Ua coordinates {W = W{xa,ya))- 

The physical role of a Coulomb-like potential is twofold. Its long-distance 
part modifies the asymptotic motion, while its short-range part strongly corre- 
lates the two-body subsystems. Merkuriev proposed to split the potentials into 
short-range and long-range parts in the three-body configuration space via a 
cut-off function ( ^ , 

Va\xa,ya) ^ Va{Xa)Ca{Xa,ya), (2) 

and 

Va\^a,ya) = Va{Xa)[l - Ca{Xa,ya)]- (3) 

The function (a is defined such that it separates the asymptotic two-body sector 
f2a from the rest of the three-body configuration space. On the region of fla 
the splitting function C,a asymptotically tends to 1 and on the complementary 
asymptotic region of the configuration space it tends to 0. Rigorously, f2a is 
defined as a part of the three-body configuration space where the condition 

IxqI < a(l -I- ly^l/a)'', with a > 0, < i/ < 1/2, (4) 
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(s) 

is satisfied. So, in i?^ the short-range part Va coincides with the original 
Coulonib-hke potential Va and in the complementary region vanishes, whereas 
the opposite holds true for Va^ . From its construction follows that if a is chosen 
big enough J2a three-body Hilbert space does not support any bound 

state This phenomena is analogous to the observation that some special 
atomic three-body systems does not have any bound states Note that for 
repulsive Coulomb interactions one can also adopt Noble's approach 0, which 
can be considered as a special case of Merkuriev's splitting, where the splitting 
is performed in the two-body configuration space. Then vi'^ coincides with the 

(s) 

whole Coulomb interaction and Va with the short range nuclear potential. 
In the Faddeev procedure we spilt the wave function into three components 

Itf') = IV'a) + m + IV'a) (5) 

by applying the asymptotic filtering operator 

\^^) = G(Hz)vi^^\^). (6) 

Here G*^'^ is the resolvent of the long-ranged Hamiltonian 

ijW =i/" + T4^ + z;«+«(')+z;« (7) 

and z is the complex energy parameter. The wave-function components satisfy 
the homogeneous Faddeev-Merkuriev integral equations 

|V'a) = GW(z)«(^)^|^,), (8) 

where Ga^ is the resolvent of the channel long-ranged Hamiltonian 

i/W^i^W+^W. (9) 

Merkuriev has proved that Egs. (^) possess compact kernels for positive E 
energies, and this property remains valid also for complex energies z — E—ir/2, 

r > 0. 



3 Solution method 



We solve these integral equations in Coulomb-Sturmian-space representation. 
The Coulomb-Sturmian (CS) functions are defined by 



{r\n) 



{n + 2l + l)\ 



1/2 



(26r)'+i exp(-6r)L^'+^(2&r) 



2/+1/ 



(10) 



with n and I being the radial and orbital angular momentum quantum numbers, 
respectively, and b is the size parameter of the basis. The CS functions {|n)} 
form a biorthonormal discrete basis in the radial two-body Hilbert space; the 
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biorthogonal partner defined by (r|n) = r~^{r\n). Since the three-body Hilbert 
space is a direct product of two-body Hilbert spaces, an appropriate basis is 
the angular momentum coupled direct product (the possible other quantum 
numbers are implicitly assumed) 

\w)a = \n)a®W)a, (n, = 0, 1 , 2, . . .) ■ (11) 
With this basis the completeness relation takes the form 

N 

1 = lim V a{ni^\ = hm 1^. (12) 

n,i/— 

Note that in the three-body Hilbert space, three equivalent bases belonging to 
fragmentation a, (3 and 7 are possible. 

In Ref. H] a novel approximation scheme has been proposed to the Faddeev- 
type integral equations 

|Va) = G«(z)lM^)^li^|^^), (13) 

(s) 

i.e. the short-range potential Va in the three-body Hilbert space is taken to 
have a separable form, viz. 



AT 

E 



nv)a v^^l i3{n'v'\, 



(14) 



where y^^^ — a{m'\va'' \n'v')fj. In Ref. the validity of the approximation were 
proved. The argumentation in Ref. relies on the square integrable property of 
the term Va'^ \''P■y)^ 7 7^ ck- Thus this approximation is justified also for complex 



energies as long as this property remains valid. In Eq. (14) the ket and bra 
states are defined for different fragmentation, depending on the environment 
of the potential operators in the equations. Now, with this approximation, the 
solution of the homogeneous Faddeev-Merkuriev equations turns into solution 
of matrix equations for the component vector i/;^ = a{nv\^a) 

where = q(?w|G'q^ InV)^. A unique solution exists if and only if 

det{[G('Hz)]"^-i;(")} = 0. (16) 

The Green's operator Gq"* is a solution of the auxiliary three-body problem 
with the Hamiltonian Ha^ . To determine it uniquely one should start again 
from Faddeev-type integral equations, which does not seem to lead any further. 
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or from the triad of Lippmann-Schwinger equations Q . The triad of Lippmann- 
Schwinger equations, although they do not possess compact kernels and thus 
they are not amenable for practical calculations, also define the solution in 
an unique way. They are, in fact, related to the adjoint representation of the 
Faddeev operator JTof . The Hamiltonian Ha \ however, has a peculiar property 
that it supports bound state only in the subsystem a, and thus it has only 
one kind of asymptotic channel, the a channel. For such a system one single 
Lippmann-Schwinger equation is sufRcient for an unique solution pl] |. 
The appropriate equation takes the form 

Gi') = G„ + G„C/"GL'), (17) 

where Gq is the resolvent channel-distorted long-range Hamiltonian, 

if, (18) 

and [/" = W + v''p + w^'-* — Ua^. The auxiliary potential Ua^ depends on the 
coordinate ya- It is defined such that it does not support any bound state and 
has the asymptotic form Ua^ ^^{ep + ^-^^IVa as j/q oo. In fact, u'a is an 
effective Coulomb interaction between the center of mass of the subsystem a 
(with charge e-y) and the third particle (with charge e^). Its role is to 

compensate the Coulomb tail of the potentials v^* -V Vj^ in f2a- If Ua^ is a 
repulsive Coulomb potential the requirement that it should not support bound 
states can can easily be fulfilled. For attractive the infinitely many bound 
states should be projected out, which leads to a non-local potential. 

It is important to realize that in this approach to get the solution only the 
matrix elements G^'-* are needed, i.e. only the representation of the Green's 
operator on a compact subset of the Hilbcrt space are required. So, although 
Eq. ( |l7| ) docs not possesses a compact kernel on the whole three-body Hilbert 
space its matrix form is effectively a compact equation on the subspace spanned 
by finite number of CS functions Thus we can perform an approximation, 
similar to Eq. (|lj), on the potential U°' in Eq. JTtI), with bases of the same 
fragmentation a applied on both sides of the operator. Now the integral equa- 
tion reduces to an analogous set of linear algebraic equation with the operators 
replaced by their matrix representations. The solution is given by 

[G«(^)]-'-[GJz)]-i-C/". (19) 

The most crucial point in this procedure is the calculation of the matrix 
elements — a{'nv\Ga\n'v')a, since the potential matrix elements v^^^p and 
C/" can always be calculated numerically by making use of the transformation 
of Jacobi coordinates. The Green's operator Gq is a resolvent of the sum of 
two commuting Hamiltonians, Ha — hx^ + hy^ , where hx^ = h^^ + Va and 

hy^ = hy^ -\- Ua\ which act in different two-body Hilbert spaces. Thus, us- 
ing the convolution theorem the three-body Green's operator Gq equates to a 
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convolution integral of two-body Green's operators, i.e. 

GJz) = ^ / dz'g^^iz - z') .g,Jz'), (20) 

where g^^iz) = (z — hx^)^^ and gy^{z) = (z — hy^)^^. The contour C should 
be taken counterclockwise around the continuous spectrum of hy^ so that gx^ 
is analytic in the domain encircled by C . 

To examine the structure of the integrand let us shift the spectrum of gx^ 
by taking z = E + le with positive e. By doing so, the two spectra become 
well separated and the spectrum of gy^ can be encircled. Next the contour C 
is deformed analytically in such a way that the upper part descends to the 
unphysical Riemann sheet of gy^ , while the lower part of C can be detoured 
away from the cut [see Fig. The contour still encircles the branch cut sin- 
gularity of gy^ , but in the e ^ limit it now avoids the singularities of gx^ ■ 
Moreover, by continuing to negative values of e, in order that we can calculate 
resonances, the branch cut and pole singularities of gx^ move onto the second 
Riemann sheet of gy^ and, at the same time, the branch cut of gy^ moves onto 
the second Riemann sheet of gx^ ■ Thus, the mathematical conditions for the 
contour integral representation of Gaiz) in Eq. ( pO| ) can be fulfilled also for 
complex energies with negative imaginary part. In this respect there is only a 
gradual difference between the bound- and resonant-state calculations. Now, 
the matrix elements can be cast in the form 

G„W ^^J^ dz'g^Jz~z')gJz'), (21) 

where the corresponding CS matrix elements of the two-body Green's operators 
in the integrand are known analytically for all complex energies , and thus 
the convolution integral can be performed also in the practice. 



4 Numerical illustration 



To illustrate the feasibility of this method we examine the convergence of the 
results for three-body resonant-state energies. For this purpose we take an Ali- 
Bodmer-type toy-model for the charged three-a system interacting via s-wave 
short-range interaction. To improve its properties we add a phenomenological 
three-body potential. Adopting Noble's splitting we have 

vi'^r) = Vi exp{-rVA2} + V2 expj-rV/Ja'} (22) 

with Vi ^ 125 MeV, V2 = -30.18 MeV, Pi = 1.53 fm, /^a = 2.85 fm, and 

^;«(r)=4eVr. (23) 

We use units such that ft^/m = 41.47 MeV, = 1.44 MeV fm. The mass of the 
a-particle is chosen as M = 3.973m, where m denotes the mass of the nucleon. 
The three body potential is taken to have Gaussian form 

W{p) = Vexp{-pyp^}, (24) 
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where = ^'^f^V = —31.935 MeV and /3 = 3.315 fm. Here stands for the 

i=l 

position vector of i-ih particle in the center of mass frame of three-a system. 
We select states with total angular momentum L = 0. In Table I we show the 
convergence of the energy of the ground-state and of the first resonant-state 
with respect to iV, the number of CS functions employed in the expansion. The 
selected resonance is the experimentally well-known sharp state which has a 
great relevance in nuclear synthesis. 

For comparison we have recalculated the results of Ref . |l^] , where a two- 
channel model has been proposed for the three-a system. Table II shows a 
good agreement between the two independent methods as far as the position 
of the resonance is considered. For the width of the resonant state we got an 
order of magnitude less which indicate an order of magnitude longer lifetime. 
The origin of the discrepancy, in our opinion, could be due to the not proper 
implementation of the asymptotic boundary conditions in Ref. [l3| . 

5 Conclusions 

In this article we have proposed a new method for calculating resonances in 
three-body Coulombic systems. The homogeneous Faddeev-Merkuriev integral 
equations were solved for complex energies. For this, being an integral equation 
approach, the no boundary conditions are needed. The procedure is an exten- 
sion of the well-established Coulomb-Sturmian separable expansion approach 
[|[ ^, In fact, the extension is nothing else but the proper analytic contin- 
uation of the contour integral representation of the Green's operator. In the 
non-Coulomb case the method is readily applicable, only the Coulomb Green's 
operators have to be replaced by the free one. 
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Figure 1. Analytic structure of g^a {z — z') Qy^ [z) as a function of z' with z = E -\-ie, 
E > 0, e > 0. The contour C encircles the continuous spectrum of hy^ . A part of it, 
which goes on the unphysical Riemann-sheet of gy^ , is drawn by broken line. 



Table 1. Convergence of the ground-state and of the first resonant-state energy (in 
MeV) of a three-a system interacting via the potential of (^2|) with increasing basis for 
the separable expansion. A*' denotes the maximum number of basis states employed 
for n and v in Eq. {UM- 
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